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Abstract. 

The existence of quasars at redshift z > 5 indicates that supermassive black 
holes were present since the very early times. If they grew by accretion, the seeds 
of mass > 10 5 Mq must have formed at 2 ~ 9. These seed black holes may 
result from the collapse and dissipation of primordial gas during the early stages 
of galaxy formation. I discuss the effects of magnetic fields on the fragmentation 
of cold gas clouds embedded into a hot diffuse phase and a virialized dark matter 
halo. The field of 10 — 4 G ejected by supernova remnants can halt cloud break-up 
at 10 4 Mq. High star formation rates keep the clouds partially ionized, making 
ambipolar diffusion inefficient. The magnetically-supported clouds collapse into 
black holes, which later spiral via dynamical friction into a central cluster with 
the total mass ~6x 10 6 Mq. As the cluster collapses, the black holes merge 
emitting gravitational radiation that should be detectable by LISA. 



1. Seed black holes in centers of galaxies 

Supermassive black holes have been inferred and detected in many galaxies, distant 
and nearby. They seem to be almost ubiquitous in large elliptical galaxies (Magorrian 
et al. 1998). The detection of quasars at redshift z > 5 (Fan et al. 1999; Stern et 
al. 2000) implies that supermassive black holes existed when the Universe was only 
tn ~ 10 9 yr old. This amount of time is not very large on galactic scales and provides 
a powerful constraint on the formation of black holes. 

An existing black hole can grow by accretion from the gaseous disk or disrupted 
stars. Accreting with a standard efficiency Em = L/Mc 2 ~ 0.1 and emitting at a 
nearly Eddington rate, Sl = L/L^dd ~ 1, the black hole doubles its mass in 
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yr. 



If the accretion is inefficient {em *C 0.1), as in the ADAF (Narayan & Yi 1994) and 
ADIOS (Blandford & Begelman 1999) models, the timescale is shorter but I consider 
the standard scenario as an example. After about ten e-folding times, the black hole 
can grow in mass by four to five orders of magnitude. Thus, starting at i# w 5 x 10 s 
yr (for a flat cosmological model with Slo = 0.4 and the Hubble constant h = 0.65) a 
10 5 Mq black hole in the high-redshift quasar progenitor can reach the mass ~ 10 9 Mq. 

However, this efficient accretion mechanism requires a seed black hole in place at 
z » 9. At that epoch large galaxies do not yet exist. In the hierarchical formation 
scenario, they assemble via infall and mergers of small clumps of baryonic gas and 
dark matter. Even though we may never know the exact route leading to black hole 
formation, we can consider various scenarios and choose the most plausible. Some 
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proposed scenarios involve stars, producing black holes either as an endresult of the 
evolution of massive stars or through a cluster of degenerate stellar remnants merging 
into a single object. These processes, relying on two-body relaxation of the stellar 
distribution, are slow and require a large number of local dynamical times, at least 
10 8 — 10 9 yr (Quinlan & Shapiro 1990; Lee 1995; Lee, these proceedings). Alternatively, 
purely gas-dynamical scenarios seem to be more efficient (e.g., Loeb & Rasio 1994). 
The collapse and fragmentation can proceed on the free-fall timescalc, < 10 6 yr at 
r = 100 pc, while the cooling time at that radius is just 5 x 10 3 yr. 

It seems that in order to form a black hole in a gas-dynamical process, two generic 
constraints must be satisfied: (1) the gas must be strongly self-gravitating; and (2) 
star formation must be inefficient, allowing a large supply of cold gas to be converted 
to a black hole. This second constraint is particularly important, since star formation 
is a competition to the black hole formation. 

In order to prevent the formation of stars, small-scale density perturbations should 
be effectively erased or damped. One way to achieve this is to keep the gas hot, with 
the thermal speed comparable to the gravitational free-fall velocity; another is to have 
an almost relativistic equation of state, where information is again transmitted faster 
than the perturbations grow. The latter criterion is satisfied in objects dominated by 
radiation pressure and by isotropic (well-tangled) magnetic field. 



2. Conservation of angular momentum 

To set the stage, consider the just virialized dark matter halo at z v ; r = 9 with the mass 
M vir w 1.8x10 10 ct? 00 [10/(l+z vir )] 3 / 2 Mq and the baryonic fraction f b = Q b /n = 0.1. 
Assume an isothermal density distribution of dark matter with a velocity dispersion 
o~dm — IOOctioo km s _1 . As the initial density perturbation grows, the gas heats 
adiabatically (T oc p 2 ^). During the mergers of pregalactic clumps, the gas is shock- 
heated to the virial temperature and follows the dark matter profile: 

For the mean molecular weight of gas with the primordial composition, ph = 0.6, 
the corresponding virial temperature is T v ; r = 3.6 x 10 5 (1 + A i t) (T ioo K, where 
/Xb(r) = - J fb(r)dr. The initial amount of angular momentum acquired from tidal 
torquing at the turn-around radius of the protogalaxy is given by the dimensionlcss 
factor, which is remarkably constant in cosmological simulations: 

J\E\ 1 / 2 
A = IJ^I = 0.05, 
GM 5 / 2 

where M is the total mass, and E is the total energy. When the gas density in the 
center becomes high enough that atomic cooling is fast (t coo \ < tdyn), the gas decouples 
from dark matter. Cooling gas contracts until it settles into the centrifugally- 
supported disk with the characteristic size Rd « Ai? v ir- 

The early epoch of a hierarchical galaxy assembly differs from the later 
times in that the infalling clumps of gas and dark matter scatter off each other, 
redistributing angular momentum and thus providing the effective "gravitational 
viscosity". Numerical simulations of mergers, with (Mihos & Hernquist 1996) 
and without (Hernquist & Mihos 1995) the effects of star formation, show that 
gravitational torques cause the gas lose a large fraction of its angular momentum, 
up to 99%, in one or two local dynamical times. The gas in these simulations develops 
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strong inflows and accumulates in the center. Any non-axisymmetric structures, such 
as a temporary bar, would speed up the process. 

3. On the possibility of a direct formation of the central object from 
gaseous collapse 

Consider now the extreme case when the gas loses most of its angular momentum in 
the subgalactic mergers. If the effective value of A is reduced by a factor of hundred, 
A e ff = A/5 x 10~ 4 , the scale- length of the gaseous disk becomes a mere 4 ctioo pc. Such 
disk is so dense that it may trap its own cooling radiation and form a supermassive star. 
This happens when the speed of photon diffusion, c/r, falls below the free-fall velocity. 
Thus, the optical depth for photon scattering needs to be at least r > 210 af 1 , A^ 2 (a 
similar scenario has been considered by Haehnelt & Rees 1993). 

There are three possible physical processes leading to photon scattering in 
the primordial gas: Thomson scattering off electrons, resonant scattering of Lya 
photons off neutral atoms, and radiation pressure due to recombination (Haehnelt 
1995). Because of the very high density of the gas, the cooling time is extremely 
short (~ 8 yr) and therefore most of the hydrogen quickly becomes neutral. The 
small ionization fraction (x < 10~ 5 ) renders inefficient both electron scattering 

— 1/2 

(r e oc x; r e ~ 0.03 cr 100 A eff ) and the radiation pressure due to recombination 
(r roc oc x 2 ; T roc ~ 0.15 (Tjqq). The dominant opacity comes from the Lya scattering 
( T Lya ~ 10 5 erfoo A^^ 2 ). The high density of the gas implies a huge star formation 
rate, 4 x 10 4 M Q yr -1 , which provides the ionizing photons that after absorption 
and recombination become the Lya photons. These photons are numerous enough to 
keep the supermassive object in radiation pressure support. (Note, however, that if 
A > 0.01 <t 10 q even Lya scattering would not be sufficient to trap the radiation.) 

The result is a Very Massive Object, which can be either a superstar or a very 
dense disk depending on the amount of angular momentum left. The structure of 
the VMO is a relativistic 7 = 4/3 polytrope dominated by radiation pressure. Such 
configuration is unstable to radial collapse due to the effects of general relativity 
( "post-newtonian instability") and implodes to a black hole, if viscous torques can 
remove the remaining angular momentum. The final mass of the black hole depends on 
how much mass the VMO can accumulate over its evolutionary timescale, ~ 3 x 10 4 yr 
(Baumgarte & Shapiro 1999). In the spherical collapse scenario, the gas accumulates 
in the center with a constant infall rate M — 3.5 a\ 00 M Q yr -1 , so the resulting black 
hole mass can be > 10 5 M Q . 

4. Effects of the magnetic field 

If the cooling gas preserves the cosmological value of the angular momentum, A = 0.05, 
thermal instability in the disk leads to the two-phase media of cold clouds in pressure 
equilibrium with hot diffuse gas (e.g., Fall & Rees 1985). The dynamics of cold clouds 
can be significantly modified in the presence of a magnetic field. A well-tangled 
magnetic field provides additional resistance to gravitational contraction and prevents 
fragmentation of clouds below a certain mass limit, the "magnetic Jeans mass" (McKee 
et al. 1993). The rms magnetic energy, $ B , would balance the gravitational pull in 
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clouds less massive than 

M. = 0.12G-V»* B « 700 (^) (A.)' M , 

where L is the size of the cloud. 

There could be a primordial magnetic field in the early galaxies, generated before 
virialization (Gncdin et al. 2000; Kulsrud et al. 1997). But additionally, after the 
first significant burst of star formation, stellar winds inject metals and magnetic field 
lines into the interstellar medium (Rees 1994). The magnetic flux of each supernova 
remnant would expand in space until the remnants overlap or fade away. For example, 
the expanding Crab nebula is ~ 1 pc in size with -B C rab ~ 4 x 10 G (Hester et al. 
1996). The number of SN remnants in the inner 100 pc would be high, Nsnr ~ 3xl0 4 , 
assuming a lifetime of the remnant of 10 5 yr. Due to the very short mean free path in 
the magnetized medium, of order 10 8 cm, the ions and electrons behave as a fluid and 
the shock fronts of the SN remnants intersect but do not penetrate each other before 
dissipating. Therefore, the magnetic flux is determined locally by the effective scale 
of expansion of a single SN remnant, or « 3 pc. The average magnetic field can then 
be as large as 4 x 10~ 5 G. 

The second generation of clouds formed within these supernova remnants would 
be magnetized. A large fraction of the material in the clouds is ionized and their 
equilibrium temperature is therefore around T « 10 4 T4 K. Such clouds can be 
approximated as isothermal and uniform. Their structure is determined by the 
pressure of the hot diffuse gas. The fraction of hot gas is determined locally as the 
amount of gas unable to radiate away the heating by gravitational infall and by other 
sources. Since in the isothermal distribution, idyn oc r, the density of the hot media 
would initially be ph(r) oc r _1 . After the onset of star formation from cold gas, heating 
by stellar winds and supernova explosions adds to the energy balance and couples the 
amounts of available hot and cold gas: 

Ph (r) = 1.5 x ID- 3 a%l r^f f^)^ g cm" 3 , 

with T = /fe,o.i /*,o.i fbfi.i A_23 Xq 05 e~ l l ta ~ 1. Here r W o is the galactic radius 
in units of 100 pc, fk is the fraction of supernova explosion energy that goes into 
the kinetic energy of the surrounding media, /* is the local fraction of cold gas that 
turns into stars on a dynamical timescale, and A = 10 _23 A_23 erg cm -3 s _1 is the 
minimum of the cooling function of the primordial gas at T ~ 10 6 K. The local 
supply of gas is consumed on a timescale to = tdyn//* with a star formation rate 
M* = 43 er 3 00 /*,(). 1 /b 1 ^aos 2 Mq yr -1 . This gives the internal cloud density 

n(r) = 1.7 x 10 3 alg r w T (^j'^ cm" 3 , 

and the Jeans mass 

Mj « 3.4 x 10 5 oZ r%t T 4 2 (I±^j M . 

If the cloud cools enough that the value of the Jeans mass falls below its mass, 
M c i, the cloud starts contracting and fragmenting into smaller, denser subclouds which 
cool further. This process can go on until M c \ ~ M$, at which point the magnetic 
pressure will resist fragmentation. Since both the gravitational and magnetic energies 
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scale with the cloud size as L~ , the cloud will remain in marginal equilibrium at 
any value of L. Thus, any increase in the external pressure will cause uniform cloud 
collapse. Unless the magnetic flux can be effectively dissipated, the inevitable outcome 
of such contraction is the formation of a compact supermassive object. 

In Galactic molecular clouds, the magnetic flux is removed by ambipolar diffusion 
of ions through neutral atoms. The diffusion timescale is 

where L is the scale of the magnetic field variation (Draine, Roberge & Dalgarno 1983). 
The efficiency of ambipolar diffusion depends on the ionization fraction, x. The cores 
of molecular clouds are shielded from external ionizing radiation (x < 10~ 7 ) and 
the magnetic field decays to allow further fragmentation and star formation. In the 
primordial clouds, the X-rays from first star clusters should be able to penetrate 
the whole cloud. The X-ray luminosity at I keV resulting from a solar mass of 
star formation per year can be empirically estimated to be <~ 3 x 10 40 erg s -1 (Oh 
2001). This leads to a significant ionization fraction, x c ~ 2 x 10~ 3 , and a long 
ambipolar diffusion time, t ad <~ 2 x 10 9 yr. Therefore, the clouds with mass M$ would 
remain magnetically supported. Also, a frozen uniform magnetic field is important for 
removing cloud's angular momentum via "magnetic braking" on roughly a dynamical 
time of the ambient medium (McKee et al. 1993). 

The dynamics of cold clouds is regulated by two competing effects, condensation 
from the hot phase and destruction in predominantly inelastic collisions. Thermal 
instability of the virialized gas leads to the formation of clouds on a cooling 
timescale. The binding energy of these clouds is relatively small; the escape velocity 
~ 2 (Ms/ Lpc) 1 / 2 km s _1 . The clouds in the disk move with a much higher velocity, 
^roi <; 10 km s _1 , and therefore cloud collisions would quickly disband them. The 
time between collisions depends on the number of clouds within the inner 100 pc. 
Assuming that the clouds develop a scale- free mass distribution dN c i/dM cx M~ 2 
with M m i n « M$ and M max « Mj, the smallest clouds would survive long enough 
for magnetic braking to operate if M$ w 3 x 10 4 r^ 3 / 4 ct^qq 3 ^ 4 M®. 

Once the magnetically-supported cloud loses enough angular momentum, external 
pressure perturbations would cause an unlimited uniform collapse to a black hole 
through an intermediate stage of a supermassive star. During the radiation-supported 
VMO phase, strong winds and outflows may reduce the final mass by roughly a half, 
although the exact value is uncertain. The VMOs close to the center would spiral 
in via dynamical friction off the gas and stars. Allowing time t<n = 3 x 10 8 yr for 
the dynamical friction with the circular velocity ~ 170 crioo km s -1 , the radius of the 
swept-up region is about r<jf = 30 pc (Binney & Tremaine 1987) and the critical cloud 
mass M$ = 1.5 x 10 4 a^J M®. This requires the magnetic field B ~ 4 x 10~ 4 G 
on the scale of the cloud, L$, ~ 2 pc. The number of M$ clouds within the central 
30 pc is ~ 780, which including the mass loss of the VMO results in the cluster of 
black holes with the total mass <~ 6 x 10 6 M®. The cluster keeps contracting until 
it dominates the gravity in the center and the dynamical friction becomes inefficient. 
The size of the black hole cluster at that point is R c i = 0.9 A0.05 pc. 

Notice the strong dependence of the black hole mass on the value of the magnetic 
field and the circular velocity of the halo, M^h oc M<j, oc B 3 a~ 7 . Small galaxies 
forming first would have a higher mass of the magnetized clouds and also require a 
weaker magnetic field to support them. 
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5. Gravitational radiation from merging black holes 

The subsequent evolution of the cluster of black holes is relatively well understood. A 
self-similar collapse occurs (Cohn 1980) until binaries form via gravitational radiation 
(Quinlan & Shapiro 1987), their eqs. (41)-(44). Because of the small number of 
objects in the cluster the initial relaxation time is short, t le i = 6.5 x 10 4 yr, but the 
initial binary capture time is long, i cap = 3.3 x 10 10 yr. Binary formation becomes 
effective when the cluster core collapses to ~ 0.02 R c i on the collisionless timescale 
t cc — 330t rc i ~ 2 x 10 7 yr. Since this time is shorter than the local Hubble time 
(t H = 5 x 10 8 yr), most of the black holes undergo hierarchical binary mergers through 
gravitational wave radiation. 

The merger of a large number of black holes would lead to a significant amount of 
gravitational radiation from the galactic center. The comoving density of halos with 
M vir = 2x 10 10 M Q at z = 9 is ~ 7 x 10~ 3 Mpc~ 3 . A coalescence of two 10 4 M black 
holes would lead to a gravitational wave amplitude h ~ 0.1 r g /Rn ~3x 10~ 20 M4, 
where r g is the gravitational radius of the black hole and Rr is the Hubble distance. 
The maximum frequency of gravitational radiation is ^ max < 0.1 c/r g ~ 1 M^ 1 Hz. 
The mergers would proceed hierarchically, doubling the mass in each event. Thus, the 
last merger of the 3 x 10 6 M Q black holes gives the maximum amplitude h ~ 10~ 17 . 

Such gravitational signal should be detectable by LISA. The early black holes 
should originate in the highest density peaks and therefore be highly clustered. They 
should appear in the same direction as the high-redshift quasars. 

I would like to thank M. Rees, M. Bcgelman, J. P. Ostriker, and E. Zweibel for critical 
discussions, and the LISA conference organizers for hospitality. 
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